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S—M' We present detailed arguments and calculations in support of our recent proposal to 

O , . identify the axion arising in the solution of the strong CP problem with the Majoron, the 

O ' (pseudo-)Goldstone boson of spontaneously broken lepton number symmetry. At low ener- 

gies, the associated f/(l)i becomes, via electroweak parity violation and neutrino mediation, 
indistinguishable from an axial Peccei-Quinn symmetry in relation to the strong interactions. 
^Nj , The axionic couplings are then fully computable in terms of known SM parameters and the 

P^ ' Majorana mass scales. The determination of these couplings involves certain three-loop dia- 

grams, with a UV finite neutrino triangle taking over the role of the usual triangle anomaly. 
A main novelty of our proposal is thus to solve the strong CP problem by a non-local term 



f-f^ ^ that produces an anomaly-like term in the IR limit. 
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1 Introduction 

In this paper we provide detailed arguments and calculations for our recent pro- 
posal to identify the axion with the Majoron [1]. This proposal rests on very 
special features of the neutrinos and their mixing which, to the best of our 
knowledge, have not been fully exhibited in the existing literature. The present 
work grew out of an earlier attempt [2] to incorporate classical conformal sym- 
metry into the Standard Model (henceforth abbreviated as 'SM'), and thereby 
resolve some of the outstanding issues of SM physics in a minimalistic fashion, 
that is, without introducing any large intermediate mass scales or unobserved 
new particles beyond the known SM spectrum other than right-chiral neutrinos 
and one new complex scalar field (and, in particular, no low energy supersym- 
metry).^ Any such proposal must in particular account for the two phenomena 
whose explanation is commonly attributed to new scales intermediate between 
the electroweak scale and the Planck scale, namely light neutrino masses and ax- 
ion couplings. In fact, it was understood already some time ago [3] that the first 
item can be achieved by taking the neutrino Yukawa couplings to be ^ 0(10"^) 
(an acceptable fine-tuning in view of the fact that the mass ratio mu/rnt is of 
the same order of magnitude). The present work, then, addresses the question 
whether one can likewise solve the strong CP problem without new mass scales 
of order > 10^° GeV. Consequently, the model that we will study in this paper 
is a minimal extension of the SM with only right-chiral neutrinos and one new 
complex scalar field (j){x), but without explicit mass terms for the fermions (and 
in particular, no explicit Majorana mass). For definiteness we will refer to it as 
the Conformal Standard Model (abbreviated as 'CSM') in the remainder. 

Apart from underlining the potential importance of conformal symmetry 
in the SM, the present proposal furnishes an (in our view) intriguing link be- 
tween strong interaction physics, where the axion is needed to solve the strong 
CP problem [15, 16, 17, 18], and the electroweak sector of the SM, where the 
Majoron arises as the Goldstone boson of spontaneously broken lepton number 
symmetry [19, 20] (a possible link between axions and neutrinos had already 
been suggested in [21, 22, 23]). This link is encapsulated in the relation [1] 

/a"-:^E-'^ (1) 

tying the axion coupling fa to the ratio of the sum of the light neutrino masses 
and a (mass)^ parameter A^^, where A^ is of order Mw or the largest of the 
heavy neutrino mass eigenvalues; a^, is the SU{2)^ gauge coupling (the estimate 
(1) is based on the formula (85), which is obtained from a three-loop calcula- 
tion and more complicated, see the final section of this paper). The weakness 
of axionic couplings can thus be naturally understood and explained, without 
adducing any large intermediate scale. 



^For alternative proposals in this direction, see [4, 5, 6, 7, 8, 9, 10, 11, 12, 13] and references 
therein. In particular, our model is similar to the so-called i^MSM first introduced in [14] and 
further studied in [7, 8], except that the extra scalar field </>(x) here is complex and not 
supposed to be identified with the inflaton. 



The main new features of our proposal can be summarized as follows: 

• If the axion is identified with the Majoron, all its couplings arise as effective 
couplings via loop corrections with the insertion of a triangle diagram 
involving neutrinos on the internal lines. There is consequently no need 
for new global symmetries beyond the global symmetries already present 
in the CSM (baryon and lepton number symmetry). 

• This triangle diagram, though superficially reminiscent of the well known 
anomalous triangle diagrams in QED and QCD, is UV finite and produces 
a non-local contribution to the effective action, which however reduces to 
the (local) anomaly-like amplitude in the low momentum limit. Hence, 
a main novelty of our proposal is to solve the strong CP problem via a 
non-local term that coincides with the usual anomaly only in the IR limit. 

• The non-vanishing result for the amplitude hinges on the spontaneous 
breaking of lepton number symmetry (as well as on electroweak symmetry 
breaking), and on the fact that parity is maximally violated in the SM. 
If lepton number symmetry were restored, the relevant diagrams would 
vanish identically even if electroweak symmetry remained broken. 

• If spontaneous symmetry breaking and the emergence of mass scales can 
be accounted for by a Coleman- Weinberg-type mechanism, starting from 
a conformal classical Lagrangian, all effects can be viewed as ultimately 
resulting from the associated conformal anomaly. 

We here concentrate on the link with strong interaction physics, while the 
discussion of photon couplings, as well as implications for cosmology, especially 
with the axion as a Dark Matter candidate, will be analyzed elsewhere. However, 
before entering into any details of the calculation it is perhaps worthwhile to 
clearly state what the result provided at the end of this paper actually means. 
We want to calculate the effective coupling of the axion to two gluons in lowest 
non-vanishing loop order. The term 'effective' in this context means that in the 
effective lagrangian we leave only massless or almost massless particles, namely 
axions, light neutrinos, gluons and photons (but the latter will not be considered 
in this paper, as we said), and therefore 'integrate out' all massive fields, whose 
masses are much larger than the external momenta, and which do not appear 
on the external legs (i.e. charged leptons and right-chiral neutrinos, quarks, 
massive gauge bosons, and all scalar fields). As we will see, there are many 
diagrams that should be considered at one and two loops, but it turns out that 
almost all of them give vanishing contribution to the coupling at hand because 
of the powerful fact that all gauge anomalies cancel in the SM. Therefore only 
the three-loop diagrams studied in the last section actually contribute to the 
effective coupling. Nevertheless, we find it convenient to divide the calculation 
into three steps associated to the different loop orders, but it is important to 
keep in mind that the results of these intermediate calculations are not the 
effective couplings of the axions to W- or Z-bosons, or to quarks, at least not 
in a proper technical sense. 



The organization of this paper is as follows. In the next section we briefly 
review the standard axion scenarios and the role of the Peccei-Quinn U{1)pq 
symmetry in them, also sketching why and how our proposal differs from the 
standard scenario. In section 3 we summarize the basic features of the CSM 
needed for our calculation; an alternative Lagrangian ('picture') giving the same 
results and based on a redefinition of the fermionic fields by a global U{1)b-l 
rotation is presented in section 4. The requisite 'neutrino technology', based 
on the systematic use of Weyl spinors and off-diagonal neutrino propagators, 
is developed in the following section, followed by a brief discussion of neutrino 
mass matrices in section 6; here we also mention some interesting physical im- 
plications, such as the possibility of light neutrino decays via emission of 'soft 
axions' (or 'soft Majorons'). Section 7 deals with the issue oi aZ mixing and the 
role of anomaly cancellations in the SM in the present proposal. The remaining 
sections address the main topic of this paper, with the calculation of the neu- 
trino triangle diagram in section 8, the axion-quark diagrams in section 9, and 
finally the axion-gluon vertex in section 10. 

A main technical novelty is our consistent use of two-component (Weyl) 
spinors for loop computations with neutrinos on the internal lines. Indeed, our 
results would be rather cumbersome to rephrase in terms of 4-componcnt (Dirac) 
spinors for the neutrinos, so this formalism affords an efficient and elegant 
method for computing higher order corrections in the presence of significant 
mixing between the neutrino components. However, for those parts of the dia- 
grams not involving neutrinos, there is not much of a technical simplification, 
so we will return to the use of 4-spinors for the quark loops in the last section. 
A short appendix explaining our notations and conventions for two-spinors is 
included at the end (see also [24]). 

2 Axions: the 'standard model', and beyond 

To contrast our proposal with current ideas and to make our presentation self- 
contained, let us first recall some key features of axion models leading to the 
standard Lagrangian governing the coupling of axions to gauge fields 

1 n/ 
>Caxion(x) = — -^a(a;)TrG^'^G^,(x) (2) 

J a 47r 

(see also [18] for a general introduction, and [25] for a review of the phenomenol- 
ogy of axions). As commonly argued [15, 16], such a coupling can be generated 
if a{x) is a Goldstone boson associated with a spontaneously broken global C/(l) 
symmetry. As we briefly explain below, this symmetry must furthermore be as- 
sumed anomalous in order to give rise to the above term. This, then, is the 
famous Peccei-Quinn symmetry U{1)pq, which is necessarily an extra global 
symmetry beyond the known global symmetries already present in the SM. The 
strong CP problem is then solved by arguing that the vacuum expectation value 
of a{x) dynamically adjusts itself to the value (a) = [17]. 



There are several axion models in the literature, all of which require extra 
and so far unobserved new particles and scales in order to render the axion 
'invisible' [26, 27, 28, 29]; see also [22, 23, 30] where the breaking of the Peccei- 
Quinn symmetry is linked to neutrino physics. Here, we briefly describe only the 
simplest model of this type. This consists in enlarging the SM by a set of heavy 
color triplet quarks Q" = [QlcnQ'ii) (where a,... are S'C/(3)c indices), which 
couple in the standard way to QCD, but do not transform under electroweak 
SU{2)u, X U{1), and hence carry no electromagnetic charge. In addition we 
introduce a complex scalar field (p which is a color and electroweak singlet, and 
couples to the new quarks via (in terms of two-component Wcyl-spinors) 

C = ^e'^^Ql^QRfia+^-c. (3) 

with the convention that a lower color index a transforms in the 3 of S't/(3)c. At 
the classical level, the full Lagrangian (including the terms we have not written) 
is invariant under the global U{\)pq symmetry 

0^e2'"0, Ql^^e-'^Ql^, QRo,a^e-'^QRo.a (4) 

In 4-spinor notation this is equivalent to the transformation Q" -^ g-i'^7 Qa^ 
whence the U{1)pq symmetry is chiral. If we now represent the scalar field as 

(f){x) = ip{x) cxp (ia(x)/v2/i) (5) 

and assume the U{1)pq to be spontaneously broken with vacuum expectation 
value (if), we can make the new quarks as heavy as we like by arranging {(p) to 
be large enough; this is, in fact, necessary in order to avoid immediate conflict 
with observation. The phase factor a(x), on the other hand, is turned into a 
Goldstone boson. The coupling (2) then effectively arises via the diagram 
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Fig.l Anomalous axion-gluon effective vertex with heavy quarks 

Observe that this diagram gives a non-zero result precisely because of the 
anomaly: if this anomaly were cancelled by further contributions from other 
fields the effect would disappear. This is the reason why in all available models 
the U{V)pQ is necessarily anomalous. 

By and large, all existing axion scenarios are variants or more elaborate 
extensions of this basic model. For instance, if one wants the axion to also 
couple to photons one must allow the new scalar field to also transform under 
electroweak symmetries, e.g. by introducing further Higgs doublets into the 
theory, as happens to be the case for supersymmetric extensions of the SM, 



thereby making room for extra global symmetries. A feature common to all 
these scenarios is that one must tune the relevant scales by hand to very large 
(or small) values in order to achieve the required invisibility of the axion. 

As a general feature we note that (2) is the only possibility to reconcile 
unbroken SU{3)c gauge invariance with the fact that a{x) is a Goldstone boson, 
which requires all its couplings to be via derivatives. This is because the gluonic 
topological density can be represented as the divergence of a local quantity, the 
Chern-Simons current JT^g, according to 



Tr G Gfj_iy — Ofj^jfjg 



(6) 



whence the coupling (2) is equivalent by partial integration to the derivative 
coupling d^aj^g. This argument and the Goldstone nature of a(x) imply in 
particular that, with unbroken gauge invariance, couplings oc aTrG^'^G^^ (or 
higher dimension couplings of this kind) can never be generated. With regard 
to the unbroken SU{i)c x U{l)em gauge sector, that is, its coupling to gluons 
and photons, the axion a{x) therefore effectively behaves like a pseudoscalar. 

In this paper we present a scenario which differs from the one above in some 
essential regards. Nevertheless the basic mechanism of identifying the axion a{x) 
with the Goldstone boson of a spontaneously broken global U{1) symmetry and of 
generating the desired coupling (2) via a loop diagram remains in force. Yet, we 
here do not assume (/'(x) to be a new scalar field; rather, we propose to identify 
it with the extra field already present in the CSM, whose non-vanishing expecta- 
tion value provides the Majorana mass term for the right-chiral neutrinos, while 
its phase is identified with the so-called Majoron, alias the Goldstone boson of 
spontaneously broken lepton number symmetry. This field is a natural ingredi- 
ent of any 'minimal' extension of the SM where this mass term is not simply put 
in by hand via an explicit Majorana mass term, and generally occurs in clas- 
sically conformal versions of the SM with right-chiral neutrinos. Consequently, 
we claim that the role of the Peccei-Quinn symmetry can be taken over by a 
known symmetry, lepton number symmetry U{\)l, whose associated Goldstone 
boson thus becomes identified with the axion. 
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Fig. 2. Contribution to axion-gluon effective vertex 



Evidently, our proposal runs counter to accepted wisdom as outlined above 
in several aspects. First of all, we here propose to cancel the strong CP term 
not by a strictly local term like (2), but by a non-local term which arises as 
a contribution to the effective action via the diagrams in Fig. 1, and which 
reduces to the local expression (2) only in the long distance (infrared) limit. In 



fact, (2) is the only term that on general grounds can survive this Umit: hke for 
the standard axion scenarios, the Goldstone nature of a{x) and the unbroken 
SU{3)c gauge invariance imply that (2) is the only possible local coupling that 
can emerge. These two principles are, however, not at all in contradiction with 
non-local interactions, which an expansion in momenta of the full effective can 
and does produce. 

Secondly, lepton number symmetry U{1)l is non-chiral, and (B — L) sym- 
metry is non- anomalous. As we will explain in much detail below, the effective 
coupling (2) nevertheless does arise in the long wave length limit due to the 
very special features of the neutrinos and their mixing. More specifically, the 
anomalous triangle involving heavy quarks (cf. Fig.l) is here replaced by a set 
of three-loop diagrams (an exemplary one of which is displayed in Fig. 2) with 
a neutrino triangle diagram at one end. As we will show, the latter, though UV 
finite, can give rise in the long distance limit to anomaly-like amplitudes due 
to the mixing of neutrino components and the fact that parity is maximally 
violated in the SM, and therefore plays the same role as the anomalous triangle 
in the usual models. One difference should be noted, though: while we make do 
here with an uncorrected quark box on the left side of the diagram, the actual 
value of the effective coupling could is affected by higher order QCD corrections, 
because as is large at small momenta (there is no analog of the Adler-Bardeen 
Theorem in the present situation). Although such corrections could change our 
prediction for fa, this is not a main concern here because the precise value of 
fa does not matter so much for the solution of the strong CP problem, as long 
as it is different from zero, as shown to be the case by our explicit calculation. 

An altogether different option starting from the present model would be to 
promote (B—L) to a ^oca/ symmetry by gauging U{1)b-l (see [31] and references 
therein, and [10] for more recent work). In this case, the Majoron would be 
swallowed by the U{1)b-l gauge boson to give it a mass, and could thus not 
become an axion. The possibility of gauging U{1)b~l is naturally suggested by 
the fact that the U(1)b-l symmetry is non-anomalous, and this is apparently 
the point of view taken in much of the pertinent literature on this subject, 
where the possibility of keeping the Majoron as an 'uneaten' Goldstone boson 
was never seriously considered, and effective vertices of the type analyzed in this 
paper were not computed (though known to exist). 

3 Basic properties of the CSM 

We refer readers to [32, 33] for the SM Lagrangian, enlarged by a family triplet 
of right-chiral neutrinos, and to [34] for an up-to-date overview. The scalar 
sector of the model consists of the SM Higgs doublet and in addition a complex 
electroweak singlet scalar field (t>{x) coupled to the right-chiral neutrinos via the 
usual Majorana term. Consequently the Yukawa couplings are given by 

+ L'e^*Y^jN^ + ^(PN^'^CY^N^ + h.c. (7) 



Here Q* and U are the left-chiral quark and lepton doublets, viz. 



Q'^\ .f h L^^\ - (8) 




while C/* and D^ arc the right-chiral up- and down-like quarks, E^ the right- 
chiral electron-like leptons, and A^* = v]^ the right-chiral neutrinos (we suppress 
all indices except the family indices i, j — 1, 2, 3; summation over double indices 
will be understood^). The Yukawa matrices 1^ ■ in (7) are general complex 3-by-3 
matrices, except for Y^^ which is complex and symmetric. $ is the usual Higgs 
doublet, and 4> is the new complex scalar field; note that this field has no direct 
couplings to the other SM fields (but such couplings will arise either through 
mixing or higher loop effects).'^ 

Since we are going to use two-spinors mostly let us spell out the relation 
with the four-dimensional spinors, in order to be completely explicit about our 
notation (see also the Appendix): 

\ uf J ' \ d'jf J ' y ef ' ' ' ^"'" ' ^ ' 

where we find it convenient to define the down-quarks with an extra rotation 
involving the CKM matrix V^^ 

d'l^^v^^^di^, 4"^y^^4" (10) 

see (11) below. Also for later convenience, we have chosen to label the inde- 
pendent neutrino Weyl spinors by different letters i^ and N, while for the other 
fermions we adopt the convention of labeling the independent two-component 
Weyl spinors by the subscripts L and R. Using the formulae of Appendix 
A, one then has the usual relations between 4-spinor and two-spinor expres- 
sions, such as for instance u*m* = u^u^^ + u^hq^uY — ^^R^\a + ^\a^^R or 

As is well known, one can use global unitary redefinitions of the fermion fields 
L^,E^,Q^, W and N^ to transform the Yukawa matrices Y^f , Y^^ and Y^^ to real 
diagonal form. To simplify the remaining (still general complex) matrices Y^^ 
and Y^ we still have three phase [= 1^(1)"^] rotations on (L%_E*), and another 
three on (Q*, [/') at our disposal, as well as the remaining unitary rotation of 
the down quarks D^. The latter can be used to represent the matrix Y^ as 

Y^ = VR^V'' (11) 

where R^ is real diagonal and V is the CKM matrix (with three angles and one 
phase); the extra factor V^ here is related to our definition (10) for the down 



^For clarity, we will, however, suspend the summation convention for flavor indices occa- 
sionally, and also use capital letters I, J, . . . for the quark flavors in sections 8fE. 

^More generally, one might replace the 'flavor neutral' complex scalar (p{x) in (7) by a fleld 
4>ij{x) transforming as a sextet under a horizontal (family) symmetry SU(3)h [22]. 



quarks. Finally, exploiting all remaining phase rotations, the neutrino Yukawa 
matrix can be brought to the form 

Y'' - Kl^P^R.K^P^ (12) 

where R^ is real diagonal, Kj^-^ are CKM-like matrices (each with three angles 
and one phases), and P,^ are diagonal phase matrices of unit determinant. The 
matrix Y^ thus represents altogether 15 free parameters (three real parameters, 
six angles and six phases). Moreover, and in contrast to the CKM matrix, Kf^'^ 
may (and K^ is expected to) exhibit strong mixing. 

The scalar sector is governed by the usual Lagrangian 

Scalar = -D^^<^^D^<^ - 9^0*9^0 - V{H, ^) (13) 

where H'^ = $^<I> and ip'^ = (jfcj). We do not further specify the potential 
here, but only make the usual assumption that the fields $ and 4> acquire non- 
vanishing vacuum expectation values by spontaneous symmetry breaking (which 
may occur either via the standard Mexican hat potential, or via a Coleman- 
Weinberg type breaking from a classically conformal Lagrangian; in the latter 
case, V{H, ip) would also contain logarithmic terms). The breaking in particular 
entails Dirac- and Majorana mass matrices for the neutrinos 

m,r-={H)Yl'^ , M,, :=(^)r,f (14) 

where H^ = $t<j). With our choice of phases M^ becomes a real diagonal 3-by- 
3 matrix. The matrix rriij stays complex, but from (11) we see that it can be 
represented in the form m = Vm^V\ with the diagonal matrix m^ = {H)R^ . 
As we already pointed out the vacuum expectation value (cp) need not be several 
orders of magnitude above the electroweak scale, but can be taken of the same 
order as (-ff); the smallness of the masses of light neutrinos is then achieved 
by taking the neutrino Yukawa couplings to be Y'^ ^ 0(10^^). ^ We also note 
that with the representation (5) and {(p) ^ the canonical normalization of the 
kinetic term is obtained with /i = {ip) and ip = {(p) + {l/\/2)ip' such that 



Id'^p'd.p'-^i 



-d^(^*d^(p ^ --d'^p'd^p' - -d'^a a^a + • • • (15) 



The field (p'{x) will not play a significant role in the remainder. 

Besides the (local) SU{3)c x 5t/(2)^ x U{1)y symmetries, the CSM La- 
grangian admits two global U{1) symmetries, lepton number symmetry U(1)l 
and baryon number symmetry U{1)b- These are, respectively, associated with 
the Nocther currents 

= e'J^'e' + D'-f''ly'-2i(f>''^^'(|) = J^ - 210^ 9^ (16) 



*In this case, the heavy neutrinos generally have masses < 0(1 TeV). Also, there are then 
no large corrections to the Higgs mass from the exchange of heavy neutrinos [35] . 



and 



J^ := ^QVO^ + ^C/Vc/^ + ^^7^^^ 



= iu*7^u* + -d'-ff^d' = -u'-f^'u' + -d''J^'d'' (17) 

o o o o 

(where by u^,d^,e^ and j^* we here denote the full Dhac 4-spinors, see above). 
Evidently both currents are purely vector-like. Furthermore, the scalar (f) carries 
two units of lepton number charge, hence lepton charge can 'leak' from the 
fermions into the scalar channel. For spontaneously broken lepton number the 
total current J^j^ remains conserved, but its fermionic part is not conserved (even 
at the classical level) because 

OXl = -i^.j {N'"Ni - N'^N^'') ^ (18) 

where the mass matrix M^ defined in (14) is non-vanishing for (ip) ^ 0. An 
important (and well known) fact is that the (B—L) current is free of anomalies, 
whereas Jj^ and Jj^ are anomalous separately. 

Finally, we write out those terms in the CSM Lagrangian relevant for the 
computation. After symmetry breaking and diagonalization of all mass matrices 
except neutrino ones, the electron and quark mass terms read 

>Cmass = - X! ("^e.e*e' + md^J'V* + m„,u*uM (19) 

i 

where, of course, {rrie-} = (we , ?7i^ , TOt- ) , and so on. Using two-spinors the mass 
terms for neutrinos and the axion interactions read 

C[2 - - U'^m'^Ni + Dlim*y^N'" + ^N'^JVPmi + iiV^(M*)*^'iV^") 
^^ (n^^M'^N^ - Nl{M*y^N^'^) (20) 



The interactions of the leptons with W and Z bosons are given by 



,(2) ^ _92 +_, _^A^ , « - -^ W/-;;l . ;=i-M"/3,A 






52 r, 
—-^t. 






COS 

while for the quarks they read 



4nt^ - -^W+V^^ul^a'^^^d%-%W;;{V^rd'l^a^''^ul^ (22) 



V2 ^ ^" "^f V2 

92 



TT^^t^ 

cos 0,1, 



2 y^La^ ^LfS "-La^ "L/3J 



with the CKM matrix V^^ . Here 52 is, of course, the SU{2)w gauge coupling. 



- sin^ t ^ , ^ 
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4 Alternative Picture 

Although we will use (7) to calculate the effective couplings, we should men- 
tion that there exists an equivalent approach (or 'picture') that emphasizes the 
fact that the phase of the complex field (/)(x) becomes a Goldstone boson after 
spontaneous symmetry breaking. Using (5) and redefining the fermionic fields 
according to 

{V{x),EHx),N\x)) ^ cxp(^-^yL\x),EHx),N\x)), 

{Q\x),W{x),D\x)) ^ cxp(^^yQ\x),U\x),DHx)) (23) 

we can replace the complex field 4>{x) by the real field (p{x). The phase a{x) 
then occurs only via its derivatives in the re-defined Lagrangian, as appropriate 
for a Goldstone boson. Defining the total {B~L) current 

J^L ■■= J^ - J^ ^ J^-L + 2 10^ d^^ (24) 

and using (5) with /i — (ip) ^ 0, this current assumes the universal form 

^^-L=J^-L-^5^« (25) 

for a current with a Goldstone boson a(x), with corresponding Lagrangian 

1 /2 

-Ccoidstone = --d^ad'^a+ -r^S^aJJs^L (26) 

Varying this Lagrangian w.r.t. a(x), the resulting equation of motion implies 
the conservation of the total current in the form 



^ Da - OXb-l = (27) 



M 

It is worth emphasizing that the redefinition (23) is also well-defined quantum 
mechanically, precisely because the (B — L) current is anomaly free in the SM, 
unlike the currents J^j^ and J^^ separately — this was our reason for including 
the quark fields in the redefinition (23). Therefore the change of variables (23) 
does not affect the fermionic functional measure, ensuring the mutual consis- 
tency of the two pictures also at the quantum level. In other words, it does not 
make any difference whether we base our calculation on the Yukawa Lagrangian 
(7) or on the vertex 9paJI^_^: the final result must be the same. 

Importantly, the conservation condition (27) for the generalized current says 
nothing about how the two contributions 9pJ^_i and Da conspire to produce 
overall current conservation by using the classical or quantum equations of mo- 
tion. All it says is that, whenever Da ^ 0, there must be a corresponding 
contribution to d^Ig_j^ ^ for (27) to be satisfied. Thus 

V2 

(^> 



Da = f-^ X ^ QiJb-l = X (28) 
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for the full classical or quantum equations of motion. At the classical level and 
with spontaneous symmetry breaking we have 

dXs-L = X = IM,, (N'^Ni ~ N^W'^) (29) 

as a consequence of (18). One easily checks that the Lagrangian (7) gives rise 
to a corresponding contribution to Da which just cancels the above term in the 
divergence of the total current to give d^Jj^j^ = 0. 

This cancellation mechanism persist at the quantum level. Here, we deter- 
mine the effective couplings of a{x) to other fields by evaluating the matrix 
elements 

(*|«5^J^-l|«)ipi (30) 

where |^) can be any (multi-particle) state involving excitations other than a, 
and where the subscript indicates that we amputate the external legs in the 
usual fashion. To get the corresponding contribution to the quantum equation 
of motion we stick on the classical fields xi, . . . , Xn associated to the particu- 
lar state 1^). Schematically, this turns the effective equation of motion for a 
(resulting from the full quantum effective action) into 

^ Da + . . . + XI • • • Xn(*|a^J^i|0)^pj + • • • - (31) 

where the new terms above and beyond (29) appear at higher orders in h and 
represent the quantum corrections to the classical equation of motion. Note that 
the terms involving x's are in non-local in general, but we are here primarily 
interested in the quasi-local approximation where we integrate out all massive 
fields, and look only at long wave-length (low momentum) excitations. Further- 
more, in the case at hand, all these corrections are due to neutrino mixing, that 
is, the non- vanishing r.h.s. of (29); if it did vanish, a{x) would be a free field. 

In this paper, we will compute various such quantum corrections to the 
classical equation of motion Da — ■ ■ ■ using the Lagrangian (7), but the above 
considerations show that the 'picture' of this section would give the same results. 
Schematically, the corrections (ordered in powers of h) are of the form 

X = X + h{WW + ZZ) + h^{ff + FF) + h^GG + ■■■ (32) 

where the letter / stands for quarks and electrons, and by F, W, Z and G, 
we schematically denote the field strengths of the associated vector bosons. 
There is nothing that forbids such local couplings to appear in the long distance 
limit. Equally important, the different terms 'kick in' at different energies, in 
accord with what we said in the introduction to this paper. For instance, at 
low momenta only the h^FF and h^GG terms are present. At higher energies 
(integrating out fewer degrees of freedom), these vertices 'dissolve' to become 
non-local, while the WW contribution is still effectively local at energies ~ Mw- 
For the aWW vertex, the equivalence of the two pictures has now been explicitly 
confirmed in [361. 
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5 Neutrino propagators 

Using two-component spinors, and after symmetry breaking, the free part of the 
neutrino Lagrangian is 






- m,,,."^Ni - mt.P^iV"-'' - -M^.N'^'N^ - -M*MN^^ (33) 



where, as before, the indices i, j — 1,2,3 label the family (the sum over which is 
understood). As already pointed out above, by making a unitary rotation on the 
fields iV*, we can bring the Majorana mass matrix My to real diagonal form, 
with strictly positive eigenvalues, viz. 



M,, 



S,,M, , 



Mj >0 



(34) 



By contrast, the Dirac mass matrix rriy remains a general complex 3-by-3 ma- 
trix. Let us emphasize, however, that the results to be presented do not depend 
on choices of phases or specific representations of these matrices, which are 
therefore mainly a matter of convenience. 

At this point we have two options. The first is to diagonalize the neutrino 
mass term (33) by re-defining the neutrino spinors 



by means of a unitary 6-by-6 matrix 



U = 



U3 Ui 



UU'^ = 1 



(35) 



(36) 



preserving the kinetic terms in (33). In this way the quadratic terms become 
diagonal, but the vertices coupling neutrinos to SM fields become off-diagonal. 
We will refer to this description as the 'propagator picture'. In the following 
section, we will discuss the neutrino mass matrices in somewhat more detail, 
and also return to eqns. (35) and (36). 

Alternatively, however, one may adopt another description (the 'vertex pic- 
ture') where the vertices remain diagonal (i.e. in the form of the Lagrangian 
(20)), but the propagators are off-diagonal. This 'vertex picture' has the ad- 
vantage that we can simply use the original SM Lagrangian for the interaction 
vertices. For this purpose we need to invert the (manifestly hermitian) operator 



K, 



i/3 











-M(5£ i 
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Q/3 



(37) 
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This operator is to be sandwiched between the niuhi-spinors [t^Q, z^", A^^, A^"] 

on the left and vp^D^ ^Nfj,N^ on the right. Because each entry here is a 
3-by-3 matrix in family space, the operator K, is thus represented as a 12-by-12 
matrix operator. It is convenient at this point to adopt the form (34), and this 
will be assumed from now on, while the matrix m is left in the general complex 
form. The matrix inversion can be performed in momentum space by iterated 
use of the formula 



A B 
C D 



A-^{A + B{D - CA-^B)-^C)A-^ -A-^B{D - CA'^BY^ 
-{D-CA-^B)-^CA-^ (D-CA-^B)-^ 

(38) 
where the sub- matrices A and D — CA ^B are assumed to be invertible (we 

can arrive at different forms of this identity by acting with the matrix 

on one or both sides of the original matrix and identifying sub-matrices in a 
different way). Identifying 



A 







, D = 



-M 



-M 



B = 



-m* 
—m 



, C = St , (39) 



we get the result for the inverse sub-matrix 



{D - CA-^B) ^ = 



V{p)*p^ M-^{p^ - m)m 



(40) 



where the 3-by-3 matrix 'Dija) is defined by 

V{p) := Up^ - m^m*)]Vr^{p^ - mV) - Mp^] = V{pf (41) 

A useful alternative form is 

+ M^/'^m^m*M-^m^mM^/'^] M^/^ (42) 

Therefore (not forgetting an extra factor of i in front) the results for the prop- 
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agator components read 



{<{^y^{y) 



3 13 I 



{vl{x)W^{y) 



K{^)Nl{y) 



(D^'^ixy^iy) 



(P'-ix^iy) 



(p-(x)7V^'^(y) 



(P'"(x)7V^(y) 



(Kixy^iy) 

(Kix^iy) 
{Kix)N^^{y) 



{Ki^Wliy) 



{N'^ix^iy) 



(N'^ixy-iy) 



{N"^{x)N'^{y) 



{N^''{x)Ni.{y) 



(2^ 



[m*P(p)*mt]'^,5fe^'P(=^-s') 



(43) 



1 + 777,* A'/ {p —m*m)'D[p)ra 



Ti^ 



p2 



(27r) 



^4 

" f r *Ti/„\*ii7 J _-ip(a;-y) 

l + mX>(p)(p2 -m^m*)M"^m''']'^ 



— -^[m*P(p)*]*^|(^^e 
(2^ 






-ip(2;-y) 



1«J ,^a/3g-ip (2:-y) 



[m2?(p)] 

[m2?(p)(p2 _ m^m*)M-^Y^ (J^e-'P^^^^) 



d'^p 

-|^ [2?(p)(p' - m^m*)M-^m^X' S^^e'^P^^^y^ 
(2tt) 

d^p 
(2^ 

d^p 
(2^ 



m^r^;^ ,e-'P(-^) 



[P(p)„. J jf,^^ 



^ [P(p)(p2 - m^m*)M-'f^^^e-^P(^-y^ 



d*p 

(2^ 

f4 



[VipYm^]-' p'^f^c 



tl*-' ^^"/^n-iP (^-2/) 



74 

-^^ [M-\p^ - m1m)V{p)m^Y' 6Je-'P^'=-y'> 
[M-\p^ - m1m)V{p)y' p'^^e-'P^'^-y^ 



d-^p 

(2^ 



^[l?(rtV]'^<5|e-'^(--^) 



where SL{2,C) indices are to be raised and lowered /ro777, the left, as explained 
in the Appendix. The bracket notation (•••) = (0|T(---)|0) is short-hand for 



15 



the time-ordered two-point function. The ie prescription for I'(p) (not written 
out here) is always such that the analytic continuation to Euclidean propagators 
works in the usual way. Note that one cannot simply use hermitian conjugation 
to check these expressions, because hermitian conjugation turns a time-ordered 
product into an anti-time-ordered product. One can check, however, that all 
expressions are consistent with the anti-commutation properties of the fermionic 
operators. 

As a crucial feature of these propagators we note the fall-off properties of the 
off-diagonal components at large momenta: unlike the usual Dirac propagator, 
these decay like |pp or even |p|'', and it is this feature which will account for 
the UV finiteness of all the diagrams that we will compute in later sections. 
Alternatively, the UV finiteness can also be seen in the 'propagator picture' 
with diagonal propagators from (35): there, the propagators have the usual 
fall-off properties, while the softened UV behavior of the diagrams is due to 
cancellations between different diagrams arising from the vertices, which are 
now off-diagonal. Of course, these cancellations, as well as the final results for 
the amplitudes, are independent of specific choices such as (34). 

Let us mention one possible application that demonstrates the utility of the 
formalism developed here, namely ncutrinolcss double /? decay, see [35] for a 
very recent discussion and bibliography. Inspection of the relevant diagram (see 
below) shows that the amplitude for this process directly 'measures' the propa- 
gator components {v^vi) and {v^J?''-) listed in (43). The |p|~'* decay of the {vl^v''a) 

and (i'^J^i) propagators for large momenta is different from the \p\~'^ behavior 
in models without right-chiral neutrinos, where the Majorana mass is induced by 
a (non-renormalizable) dimension-5 operator ~ e^'^Xy ($-^eL^)($^eLo) -l-h.c. 
In the neutrinoless double beta decay the external momenta are very small, but 
for larger momenta this behavior could be used in principle to discriminate be- 
tween our model, and one where the left-chiral neutrino is treated as a Majorana 
particle. 



(A,Z + 2) 



{A,Z) 




■ e 
Fig. 3 Neutrinoless double beta decay 



6 Light vs. heavy neutrinos 

Although our main result docs not depend on such choices, it is occasionally 
useful to represent the fermion neutrino matrices in a specific form; we refer 
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readers to [37, 38] for up-to-date discussions of neutrino masses and mixing. 
In particular, assuming the real diagonal form (34) it is not difficult from our 
formulae to recover the usual seesaw mass formula [39, 40, 41, 42] from the 
poles of the propagator matrix 'D{p). More precisely, the mass eigenvalues are 
obtained by replacing p^ in (41) with the parameter A, and demanding 



det {X-m^m*)M-^{\-m'im) - MX 







(44) 



In this form, the small and large eigenvalues (whose reality follows from the 
manifest hermiticity of the mass term in (37)) are still 'entangled'. For this 
reason we factorize the total mass matrix by ordering it in powers of the 'small' 
matrix m as follows,^ 



det 



(XM^^ - M - m^m*M~^ - M~^m'^m + M" 



''^■m^m*M~ 



nvra 



•(A 



M ^m^m*M '^m)m + ' 



= 



(45) 



where dots stand for higher powers of the small matrix mM^^ . It is then obvious 
that, in very good approximation, the determinant equation factorizes into a 
product of a factor det(A — M^) yielding the large eigenvalues (Mf , M2,M^) in 
(34) for the heavy neutrinos, and a second factor for the light neutrinos. The 
latter can be cast in the form 



dct(A — m'm) = , m := mM' 



-1 T 

m 



= m 



(46) 



yielding the mass eigenvalues of the light neutrinos (here we have used the 
fact that, for arbitrary square matrices A and B, the matrices AB and BA 
have the same eigenvalues). The (mass)^ matrix mm^ thus represents a matrix 
generalization of the usual seesaw formula. 

Instead of looking for the poles of 'D{p), we can also arrive at this result by 
direct consideration of the neutrino mass matrix. The parametrization of mass 
vs. interaction eigenstates of the light neutrinos is usually given as 



(47) 



where U is a. unitary matrix. In our case this formula is incomplete, due to 
the admixture of the heavy neutrinos. To spell out the precise relation between 
the 3-by-3 matrix U and the unitary 6-by-6 matrix U introduced in (36), we 
substitute the redefined fields from (35) into (33), demanding the mass term 
[= second line of (35)] to be diagonal in the redefined neutrino fields. Keeping 
in mind that this redefinition docs not mix spinors with dotted and undotted 
SL{2, C) indices, we obtain the condition 



l^e \ 




VX 


Vi, 


= U 


V2 


Vr 1 




\ v-i 



u* 







m 
M 



wt = 



m 






M' 



= U 







ra 
M 



w 



(48) 



^We here assume (mainly for simplicity) that all eigenvalues of M^ arc large in comparison 
with those of rri^m; for other mass patterns these arguments may have to be revisited. 
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where m' and M' are real diagonal 3-by-3 matrices. This implies 

,'2 



m'" 
M'2 



= U 



U 



' m*' 
m) M 


U^W 


■ 


771 

M 


m*w7' m*M 


U^ 


Mrri^ mtm + M^ 





U^ 



(49) 



To relate U to U, wc make the ansatz 
U 






U 

Um 



n 



(50) 



where U and Um are both unitary 3-by-3 matrices. Inserting this ansatz into 
(48) and expanding in powers of the 'small' matrix mM~^ up to second order, 
we obtain 



n = 



1 - ^m*M~^m'^ 



M- 



-1 T 

m 



\M-^m'^m*M- 



+ [{mM-^y 



(51) 



and, up to diagonal phase redefinitions, the matrices U and IAm are determined 
by the conditions 



m 
M' 



-WmM-^m^U'^ 



U 



M 



1 



M +-wF'm*M-^ 



1 



M-'^Tv)ra]ulj (52) 



Multiplying these matrices by the complex conjugate matrices and using the 
reality of m' and M', we see again that U diagonalizes the generalized seesaw 
(mass)^ matrix m^m. In other words, we have now rather explicit expressions 
for the eigenvalues of mass matrices of light (m') and heavy (M') neutrinos, as 
well as for the unitary matrix relating the 'propagator picture' and the 'vertex 
picture'. In lowest approximation the matrix lA in (47) is equal to Ui in (36). 

Obviously, the mixing of light and heavy neutrinos has several interesting 
physical implications. For instance, the heaviest among the light neutrinos can 
decay into the lightest neutrino via the emission of 'soft axions'. This process is 
possible because in our scenario the axion is expected to be almost massless [for 
instance, with a mass Waxion = ©(lO^^eV)], and thus much lighter than even 
the lightest neutrino. 

7 aZ mixing 

We are here concerned with the effective couplings of the 'invisible axion' a{x) 
to 'visible' SM fields. Such couplings must arise through loop diagrams (with the 
exception of couplings to neutrinos through the mixing described in the previous 
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section) , as there are no direct couplings at the tree level, a fact which according 
to our proposal can explain the extreme smallness of the axion couplings to 
standard matter. Assuming {ip) ^ we now fix the free normalization parameter 
in (5) once and for all to the value ji — (if), as in (15), in order to obtain the 
canonical normalization for the kinetic term of a{x) in the classical lagrangian. 

However, before we proceed to the actual computation we need to discuss 
the mixing of a with Z bosons. The result described in this section shows that 
in principle, there can be large effects (in this case at one loop). As we will see, 
these fail to contribute to the axion-gluon coupling only because of the vanishing 
chiral anomaly of the SM, but they can nevertheless dominate in other processes. 
There is a similar mixing between the Z boson and the scalar f'{x), as well as 
a mixing between a{x) and the standard Higgs boson, but these couplings turn 
out to be suppressed by an extra factor of the light neutrino masses, and can 
thus be neglected. 

The mixing of axions with gauge bosons, which may arise at one loop or 
higher loop orders, may in principle occur between a{x) and any neutral gauge 
field Af^{x) at the quadratic level, leading to extra terms ex d^aA'^ in the effec- 
tive Lagrangian. However, gauge invariance immediately forbids such couplings 
if the gauge symmetry is unbroken, whence a{x) cannot couple in this way to ei- 
ther photons or gluons. For the broken sector, this argument does not hold, and 
by charge conservation, we are therefore left with possible quadratic couplings 
of a to the Z boson, 

-Cmix = e^^aZ^' (53) 

where the parameter e is of dimension one. 

We now determine the mixing coefficient e at one loop (with right-chiral 
neutrinos in the loop), and show that e is proportional to the sum of the light 
neutrino masses. The relevant diagram is shown below, and involves the off- 
diagonal (Niy) and (Ni?) {{Ni') and (Np) in the second amplitude) components 
of the neutrino propagators in the loop. 

k + q 



7 -^ 


k 


a 


Fig. 4 One loop 


axion- 


-Z^ mixing 



This diagram gives rise to the following Feynman integral which can be 
evaluated in the standard fashion; note that, while naive power counting would 
suggest the presence of divergences, this diagram is actually UV finite because 
of the fast decay of the propagator components {Niy) , (ND) , {Ni') and (ND) at 
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large momenta. 



■ 92 -^Ai/Ji 






• (i[M-i((/c + qf - m)m)V{k + q)m^Y'5%l)] ] 



(54) 



.92 



■E 



d^fc 



Tr (ct^ct"^) X 



2V2(^)cos6'vF ^ J (27r)4 
X m*X>(fc)*[(fc^ - m)m){k + g)^ - fc,^((A: + qY - m)m)]'D{k + q)m 
d^k 



92 



V2{ip)cos0wJ (27I-) 



rTr 



.92 



V2{(p) COS ^IV 



E( 



TO TO, 



TO*X'(fc)*[g^(P - mV) - fc^(2fcg + q^)]Vik + q)m^ 
d*k [kf'{2kq + q^)-qf'k^]Mf 



(27r)4 /c2(fc2 _ M2)(fc + g)2[(fc + g)2 _ M2 



We are here interested only in the result for small axion momentum q^ , in which 
case the integral can be worked out to be 



AM^,,{q) « ig. 



32 



= -9' 



\/2((p) COS^TV 
p .92 



^(m'^TO*)" / 



d^fc (2fc^fc'^ - r]f'''k^)Mf 



(27r)4 [fc2(fc2~M2)]2 



(47r)22y2(93)cos0VF 



El 



,y|2 



(55) 



Remembering (14) and making the (reasonable) assumption Ym ^ 0(1) we see 
that 



,ij\2 



E l'"- I V^ 



^,^3\2 



M 



E 



my 



(56) 



whence we arrive at the claimed result e ^ "^niy, that is, the mixing is pro- 
portional to the sum of the light neutrino masses. For the mixing between Z^ 
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and the scalar ip' the two contributions to the above integral would appear with 
opposite signs, leading to a further cancellation, with a mixing parameter of 
order ©(m^) (which hence can be ignored). 

The aZ mixing described above can lead, via subsequent Z couplings, to 
further couplings of the axion to other SM fields like quarks and leptons. When 
considering physical effects such as axion cooling in stars from axion emission 
from leptons or quarks, these contributions must be taken into account. For 
the same reason, one would also expect the above mixing to contribute to the 
effective coupling of axions to gluons via the anomaly diagram 



a 




9 

Fig. 5 Potential two-loop contribution to axion-gluon vertex 

which, if non- vanishing, would be more important than the three-loop compu- 
tation we are going to perform in the final section of this paper! However, here 
we meet a pleasant surprise, which will simplify our subsequent considerations 
substantially: when summing over the quarks in this diagram one obtains a van- 
ishing result because of the cancellation of all gauge anomalies in the SM. In 
other words, when determining the effective coupling of the axion to gluons we 
can ignore the above mixing (and similar diagrams). The same conclusion holds 
for the couplings of the axion to photons, when all SM fcrmions arc summed 
over. The important fact is therefore that the non- vanishing effective couplings 
arise solely via the three-loop diagrams with an attached neutrino triangle which 
we will work out in the remaining sections. 

8 The aW^W triangle at one loop 

Our main proposal relies essentially on a new effect producing an 'anomaly-like' 
amplitude from a triangle diagram involving neutrino triangles, which in turn 
gives rise to the effective coupling of a{x) to VT-bosons; the relevant diagram 
is depicted in the figure below. The similarity of this triangle diagram with the 
well known one producing the triangle anomaly is obvious, yet the anomaly-like 
effect here is not due to a linear UV divergence, but rather to the mixing of 
the neutrino components. Technically speaking, one reason for this is that the 
neutrino propagators involve cr-matrices rather than 7-matrices, and the trace 

^M^pA ^ rj.^-^^^-p^x ^ 2{r]'"'f]P^ - ry^^ry'^^ + tj^^tj^'p + ie^"^"^) (57) 

over cr-matrices generates both parity even as well as parity odd terms, unlike 
the corresponding 7-matrix trace Tr'f^'f'^jPj^. 

In this section we take the first step towards our final goal of determining 
the effective axion-gluon coupling by calculating the aW^W^ vertex; this result 
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will then be used as an input in the calculation of the axion-quark diagram in 
the following section, which in turn will be substituted in the final step into 
the three-loop diagram yielding the effective axion-gluon vertex. Once more, we 
emphasize that all diagrams in this and the following sections are UV finite, 
even though naive power counting might suggest otherwise. Let us also point 
out that there are similar diagrams with two external Z-bosons, where the 
triangle is 'purely neutrino' (that is, all internal lines are neutrino propagators). 
In accordance with our basic strategy for computing effective couplings outlined 
in the introduction, we will however disregard these diagrams, because they 
necessarily contain a light neutrino propagator {{w) or {vp) or {iyv}) on the 
internal line connecting the two Z-boson vertices. 



P+q 




Fig. 6 Axion-W^ -W vertex 
The above Feynman diagram corresponds to the integral 



-i-^aw(P'9) 






(eUe-l<,^(fc-rt)(-i^. 



_._32 _,.A,,32 I X 
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)l3ia2 



•(i(^+^)„3AjP(fc + (7)m^P) + 

• (i[M-i((/c + g)2 - w)m)V{k + g)™^]^'**?;^;')] } = 
d^fc {k-p)x 



■ 92_ ua-2l32 



(58) 



ffi 



■E 



('^^)"^''V^^^,(a'^) 



l^A"2&^P 



272(^)^7 (2^)4(fc-p)2 
TO*2?(fc)*[(fc — 'm^'m){k + q)p — kp{{k + q) — 'nrin)]'D{k + q)in 
d^fc (fc - p) 



^2"! 



5i 



E 



-Trial" a^a-'aP} 



2^/2(^)^y (2^)4(fc-p)2-m2_ 
m*V{k)*[kp{2kq + q^) - qp{k^ - m^m*)]X>(fc + g)™^ 



(59) 



As before, we are here interested in the leading terms for small axion momentum 
qi" and small m, in a basis where M is real diagonal; in this approximation the 
integral simplifies to 



^^aWwiP^^i) ~ V 



9l 



2V2(^) 
d^fc 



E 



m*J? X 



Trja'^a'^aV} 



(fc-p)^ {2kpkr-riprk^)M^ 



(60) 



(27r)4 (fc - p)2 - m2^ [fc2(fc2 ~ M2)]2 

This integral can be evaluated by means of Feynman parameters, with the result 



-'^^aWwiP^l) 



9l 



327r2V2(^) 



Vlm^^'pMjt^-^'^'' / dy /" dxx{l-x)i 
Jo Jo 



[(-1 + 2y)pxrirp + (1 - y){PTVp\ +PpVtx)] 
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-y(l — y)p2 + (f — j/)m2 + yxM. 

2y(l - yfpxj-p^Tirp + 2prPp) 
[— y(l — y)p2 + (1 — y)m1. + yxM'^Y 



(61) 



where we made use of the definition (57) (note that the denominator in the 
integrand will become positive definite after Wick rotation to Euclidean mo- 
menta pi"). Although we do need the full expression (for large p^) below, it is 
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nevertheless instructive to specialize this result to small p to get 

„2 

^^2 J^lm'^^Mft^'^^Pq^x (62) 



^ f\ f\ a:(l - x)y^[{-l + 2y)pxVrp + (1 - y){PTVpX +PpVr\)] 
Jo Vo [{l-y)ml+yxMn^ 



a2 



SirV^if) 



ft^^i^PqT J2 



m 



ii|2 



Mi 



Mf{Mf - 5m2j log ^ + (M| + 3ml){Mf - m^J 

PA r/rp 



-(PrVpX+PpVrx)- 



6(M2-to2 )3 



Mf{Mf + ml) log ^ - 2Mf{Mf - m^J 



6(M|-m2.) 



The singularity in this expression for rrig^ = M^ is spurious. 

Let us pause to put this result in perspective. The trace over a-matrices 
contains both parity even as well as parity odd terms. The former lead to non- 
gauge invariant contributions for the amplitude, proportional to q'^p'^ + q'^p'^ 
and (q''Pp)ri^'^ , respectively. As we already pointed out, such contributions are 
to be expected because the electroweak symmetry SU{2)w x U{1)y is broken. In 
addition we get a gauge invariant anomaly-like term oc £'^^'"^Ppqo ■ The presence 
of both these terms is in accord with the fact that parity is maximally violated 
in the SM. By contrast, for gauge bosons associated with an unbroken gauge 
invariance on the external lines, only the gauge invariant anomaly-like contri- 
bution can survive by the general arguments given in section 2. Consequently, 
insertion of the above triangle as a subdiagram into a higher loop diagram with 
external photons or gluons will yield only the gauge invariant anomaly-like am- 
plitude, as we shall explicitly verify. In particular, for the axion-gluon amplitude 
we can anticipate that the result will be proportional to ie^^^'^qxPa, as SU{3c) 
remains unbroken. This is the core effect which justifies our claim as to the 
emergence of the effective coupling (2). In the following sections we will verify 
this claim by explicit computation. 

9 Axion-quark diagrams at two loops 

The next step is the the calculation of the axion-quark diagram, which is given 
by the following two-loop diagrams: 
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p + q 




p + q 



Fig. 7 Axion-up and axion-down quark vertices 

Since this diagram involves the gauge boson propagators on the internal lines 
we have to specify the gauge. In the so called R^ gauge the propagator reads 



{W+^W-''){k) 



1 



M^ 



V""" + (C - 1)- 



kf'k" 



ml 



(63) 



In the full calculation one must also include the diagrams with charged Gold- 
stone bosons from the Higgs doublet whose propagator behaves as 1/(A;^— ^M^), 
and whose contribution vanishes only in the limit ^ — > oo. The calculation then 
gets very involved even by comparison with the formulae of this and subsequent 
sections (the full calculation and all details will be given in [43]). The estimate 
shows that the contribution of these terms is of the same order of magnitude 
as what we calculate because of the large Yukawa coupling of the top quark. 
Therefore the final result that we get by neglecting both k^k" parts of the prop- 
agators and the contributions from charged Higgs particles is only an estimate 
of the value of the actual result. 

To make the formulae more transparent, we will use capital indices I, J, . . . 
for the quark flavors in the remainder, and now also write out the explicit sums 
over them. With this convention, the above diagrams correspond to the following 
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Feynman integrals 



K 

X 



{2n) 



-i^lv^K 



'V2' 



i(y - H2) 



ih + qf - M^ kl - Ml 



w 



■^aJ^ylKryUKJ f ^ 

2 V ' (2^^)' 



ip-k2f-ml^^ V '\/2 
ip - k2)T 



_;EL(y^)KJ^^ 



(27r)4 {p - k^Y - m\^ {k2 + qf - M^, kl - M^ 



w 



and 

-'^Miid{p,q) = 



(64) 



E 

K 

X 



.92 



_jjy2(0)/K-^ 



(2^)4 V ^^'^ ' V {k2~pY-ml^ 



-Ah-i) 



{k2+qr-Ml,kl-Ml, 
'2 r d*k 



i|E(^^r^^V^ 



'2 (^2 - P)t 



V2 



1 



1 



^ . ,_,r)4(p-fc2)2-m2^^(fc2 + 9)2-M^fc|-M2, 

xa^a^a''x(-iM^^^(fc2,q)) (65) 

Using the ■M'^'(y]y result already calculated before together with 

CT^^a^ X t^^^P = 86^aP (66) 

we get, as always for small g^, 

-^Mii^{p,q) = 



-iqPd" 



9\ 



Jo Jo 



d^fc ~{k-p)^ 1 

(27r)4(fc-p)2-m2^(F-M2,)2 ^ 

[(-1 + 2y)kxr]p„ + (1 - y){kpVaX + fc^^ypA)] 
[-y(l - y)A:2 + (1 - y)m2. + yxM'jY 

2y{l-yfkx{-Prjp„ + 2kpka) 
[-y(l - y)k'^ + (1 - ?/)»"4 + V^^lY 



(67) 
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and 



-iqPd^ 



9\ 



d^k {k ~ p) 






(68) 



y^ X 



^2 - Af 2^)2 



(27r)4(fc-p)2-m2^(fc2 

[(-1 + 2y)kxr]pa + (1 ~ y){kpT]a\ + kgijpx)] 
[-2/(1 - y)F + (1 - y)TO2^ + ya;M2]2 

2j/(l-J/)2fc;,(-fc2^p^+2fcpfc,) 



+ 



-y{l - y)k'^ + (1 - y)»7i2^ + yxM'jf 
Again we employ Feynman parameters to obtain 

^ qP^- 52 y T/^^(T/t)^'V-'f M2 /"'dx /"'dy /"' dz C 

64y27r4(^) ,.4- ^ ^ ' ' ^io io io Jo 



f^{x,y,z,t)r]p„ 



dt 



/^(x, y, z, t)p^f]pc, + /g (x, j/, z, t)ppPa 



/^(x, y, z, t)p'^rip„ + f^{x, y, z, t)p^PpPa 



+ 



12.^1 



x,yiZ,t;p)Y 



(69) 



and 



i-^add(P, 9) 



dz / dt 



-a 12 Y,yUIKyKJ^^J^2^2 f dx f dy f dz f 

ff{x,y,z,t)'i]p„ f^{x,y,z,t)p'^T]pa + fi{x,y,z,t)pppa 



ijk{x, y, z, i; p)]2 [M2ij-K (x, y, z, i;p)]3 

/|(a:, 2/, Z, t)p^Vpa + /5^(x, y, Z, t)p^PpPa \ 
[M2,,K(x,y,z,i;p)]4 J 

where various functions depending on the Feynman parameters are defined by 

f^{x,y,z,t) = ^x(l-x)y^z(l-z)t^(-l + 3y + 3tz) 

f^{x,y,z,t) = x(l-x)y4(l-y)z(l-z)t4(l-t)[-l + 2y + tz(l-y)(-3 + 5t)] 

/3"(x,y,z,t) - 2x(l-x)/(l-y)2z(l-z)t4(l-f)[l + fz(3-4t)] 

/4"(x,2;,z,i) = -3x(l-x)y'^(l-y)3z2(l-z)f^(l-t)3 

f^{x,y,z,t) = Qx{l-x)y'>{l-yfz\l-z)t'>{l-tf (71) 
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for the up-like quarks, and 
/f(x,2/,z,t) - -^x{l-x)y^z{l-z)t'['^-?,y-2,{l-y)zt] 
f^{x,y,z,t) = x{l~x)y^il~yfz{l-z)t\l-t)[l-^z{l + t)] 

f^{x,y,z,t) = x{l^x)y^il^y)zil-z)tHl-t)[y + 2il-y)zi^-2t)] 

f^{x,y,z,t) = 

f^{x,y,z,t) = 3x(l-x)y5(l-y)V(l-z)t5(l-t)3 

for the down-Hkc quarks. We also introduced the shorthand notation 

M^Kix, y,z,t;p) := xyztMf + (1 - y)ztml^ + y{l - y){l - z)tM^ 
+y(l - y)(l - t)ml^ - y{l - y)i(l - t)p^ 

and 

M\jKix,y,z,t;p) := xyztM^ + {1 - y)ztml + y{l - y){l - z)tM^ 

+2/(1 - y){l - t)ml^ - y{l - y)t{l - t)p^ 

Notice the difference between up-like and down-like quarks in these expressions 
(apart from the different masses of up- and down- like quarks): although the 
integrals (67) and (68) look almost the same, the indices on q^, a^ and the loop 
momentum k^ are contracted differently 
For small p we arrive at 

-iMiL{p,q)^ (72) 

x{l - x)y^z{l - z)t^{-l + 2,y + 3tz) 
"" [xyztM] + (1 - y)ztml^ + y{l - y){l - z)tM^ + y(l - y){l - t)mlj^ 

and 

-\MiUp,<l)= (73) 

x{l - x)y^z{l - z)i3[4 -iy- 3(1 - y)zt\ 
"" [xyztMf + (1 - y)ztml + y{l - y){l - z)tM^ + y{l - y){l - t)mlj^ 

In good approximation we can now put mup — rndown = (there are no IR 
divergences); then the sum over K can be performed, and by the unitarity of 
CKM matrix the amplitudes become proportional to S^"^ , i.e. flavor diagonal. 
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The non-degeneracy of the quark masses makes possible quark flavor change in 
the interaction with the axion, but the off-diagonal amplitudes are suppressed by 
factors of order 0{m'^^^j.^^/Mj). Note that the difference between the integrals 
involving up and down quarks, respectively, is due not only to the CKM matrix 
and the different quark masses, but also to the different isospin and the different 
topology of the diagrams; this leads to different formulae for the two cases. Also, 
it appears that both the gauge invariant and non-invariant parts of the aWW 
amplitude are important, as both of them contribute to this amplitude, as we 
can check in relations (66). 

10 Axion-gluon vertex 

After these preparations, we are ready at last to tackle the final part of the 
computation, which will yield the coupling of a{x) to gluons. In leading order, 
this coupling is given by the set of three-loop diagrams depicted below. The 
first set consists of the following diagrams with insertions of the axion-quarks 
diagram determined before, with either up-quarks running in the loop 



^fc-l-g 





Fig. 8 Axion-qluon vertex with up quarks in the loop 
or with down-quarks: 





Fig. 9 Axion-gluon vertex with down quarks in the loop 
However, there are also the following non-planar diagrams: 
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9p, 'TSTTBTSTSTS^ 




^ "^ ki + k2+p 

-ki - q 




9^ nrsTTBTTir^ 
l)%q 



""-k^-^-p-q^ 



Fig. 10 Axion-gluon vertex with nonplanar loop 

We first of all see that the part of the total (summed) amplitude linear in q 
is antisymmetric under the simultaneous exchange fi ^ v, p ^ —p. This means 
that the only possible tensor structures in the amplitude are either proportional 
to p'^q'^ —p'^q'^ or to e^'^^^pxqa- As we already explained, the first structure can 
be excluded by gauge invariance, which leaves only the second contribution. Of 
course, this claim is confirmed by the explicit calculation. Hence, the effective 
interaction of axion and gluons for sufficiently small q is indeed of the form (2), 
with a non-vanishing (but small) coefficient. 

The full integrals are now very cumbersome, so we try to present the result in 
a compact form. As already mentioned in the introduction, we revert to 4-spinor 
notation for the loops not involving neutrino lines. For the first two diagrams 
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we have, to leading order in q, 
— iA^? ^"n (through up quarks) = 



(agg) ^ 



-iS^" ^ 



E 






Tr 



If, + ^ + m 



[k + q] 



"^ -yA" 



9 9 



{k-pf 



uj ^ If + m 



m 



ui 



fc2 



mt 



■[-^MiUk,q 



+ (fi ^ ly , p ^ -p - q) 

_■ pgab 44 y- |y/K|2|^.,|2^2x 



X / dx Ay I dz dt 
Jo Jo Jo Jo 



Tr 



m 



MJ_ U 



^-. 






k^ — m'^ {k — p)^ — m^ k^ — rn^ 

fi{x,y,z,t)-qp„ f^jx, y, z, t)k'^i]p„ + f^{x, y, z, t)kpk, 



llj,{x,y,z,t;k)]^ [M2^.^(x,y,^,t;fc)]3 

f"{x,y,z,t)k'^T]pa + f^{x,y,z,t)k'^kpk„ 



i^{x,y,z,t;k)Y 



-P) 



(74) 



where the expression for Mi'lu{k,q) obtained in the foregoing section has been 
used. The remaining momentum space integral can be done in the standard way, 
introducing yet another Feynman parameter u (so there are now altogether five 
Feynman parameters x,y, z,t,u). Packaging the previous results into various 
new functions we get, after some calculations. 



-iAlf^'^s (through up quarks) 



i£^'^^>4g.^° \, £fi^, , E iV'^flm'^^Mfx 



32V2(27r)6(v3) 



(75) 



X I dx dy I dz I dt j 
10 Jo Jo Jo Jo 



Fi{x,y,z,t,u) 2 



where 



-^W = xyztuMf + {l-y)ztuml^+y{l-y){l-z)tuM^ 



w 



+2/(1 - y)(l - t)uml^ + y(l - y)t{l - t)(l - u)ml^ (76) 
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and the function 
Fi{x,y,z,t,u) = u{l-u) 



(-4u)/r + [(-2 + 6u)f^ - m 
+ [(4 - 8u)f^ fi 



y{l-y)t{l-t) 



[y{l-y)t(l-tW 



(77) 



itself depends on the previous Feynman parameter functions (71) and (72), but 
now with an extra dependence on the fifth Feynman parameter u. 
Similarly, the third and fourth diagrams give 



-iAl?^'^N (through down quarks) 



ie^'-^^^g.-J"" 



9 4 

.93.95 



32^2(2^)6 (^)^f7-^ 



J2 iV^^Plm^JfM^x 



(78) 



X ;'dx f'dy !\z l\t /'d.:^M^#^ + 0(p2) 











JQ 







[Mj -12 



'ijKI\ 



Here we have subsumed the previous results into the functions 

MJjKi = xyztuMj +{l-y)ztuml^+y{l-y){l- z)tuM^ 



w 



+y (1 - y) (1 - t)u ml+yil- y)t{l - t) (1 - u)ml^ (79) 



and 

F2{x,y, z,t,u) = u{l — u) 



(-4u)f^+[{~2 + 6u)f^-f^] — 



y{l-y)t{l-t) 



^^^'-'-^f'-f'^Mi^l 



[y{l-y)t{l-t)r 



(80) 



Finally, the last two diagrams must be computed directly. To leading order in 
q, they are given by formulae 



-iAl? ^"^N (nonplanar diagrams) 



= -6' 



2 2 
ab9z92 






d^ki f d^k2 



4 / (27r)4 

1 1 



(-i-^aw(^2,(Z)) 



A-2 - M2 



(fcl +p)2 - TO^^^ ;,2 _ j^2^^ (^^^ ^ ^^■)2 _ ^2^ (^^ _|_ l,^ _|_p-)2 _ ^2^ 

xTr|(^i+y + m<i^)7^(^i+mrf^)7«PL x 

X (^1 + ^2 + m„j7''(^i + ^2 + y + mu.hxPL 
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^ 4 -^l ' y (271)4; (2^)4 I '^aWwi l^2,q)) 



111 1 



2-m2^ 



(fci +p)2 - m^^ fc2 - mj^ (fci + fc2)^ - m2^ (fci + k2 + p) 
xTri (-^1 - :^ + md^-)7''(-^i + ?™d^)7K-PL x 

x(-^i -h+ TO„j7^(-^i -h-^ + mu.hxPL [ = 

.2^2 /■ ^41,. r rlil^^ f I \ 2 



,ah9m^\,rIK\2 f d ki f (£k2 



,j, . c^^rJ (2^)' \ki-My^ 



)4 y 1,2 _ /)^2 
1 1 1 



(fci + p)2 - m^^^ fc2 - m^^ (fci + fc2)2 - to2^ (fci + fc2 + p)2 - m^^ 

X T,^;:(fci,fc2,p) [-iX2^,^(fc2,(z)] (81) 

where 

Tlf;;{h,k2,p) :- (fci+p)"(fci+fc2+p)''x (82) 

X [Tr{7„7^-^i7«(-^i + ^h'^lplxPL} - (m ^ '^)] 
+ m2^(fc+p)"[Tr{7o7^^i7K7"^7APL} - (m ^ /^)] +0(m2^) 
After a tedious calculation we obtain the result 
—iAlJijI'^'' (non-planar diagrams) ~ 

2 4 

x/"d./'dj,rd3/"dt/"d„-'^*'»-^-''"' 



JO JO JO JO JO [MIki? 

where 

Fsix, y, z, i, u) == 3a;(l - x)y^z{l - z)t'^u'^{l - u) {2t'^u{9 - 4m)(1 - y)z 
-t[z{l - y){-8u^ + 15u + 1) + 2(1 - u)(l - 2y)] - u{l -y){z + 2)} 
In total, we thus arrive at the final result equivalent to (2) 

-iX^^^-^ (total) = -[-Aj-e^-P-p^q,6'^'> (84) 



where the axion coupling is given by 

128V27r4(^) ^f-.' 111^ 



92 \' It/JA|2|_I7|2a^2 ,, /gg-v 



' [' [' f' f' ( Fi -F2 + -F3 ' 

X / dx / Ay I dz At I Au\ -rr-p^ tt H 

/o Jo ^ Jo Jo Jo U-^k/]' [M^KI?, 
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This integral cannot be evaluated in closed form, but we can easily get a nu- 
merical estimate. First, to recover the order of magnitude estimate (1) quoted 
in the introduction, we take the heavy neutrino masses degenerate, that is, 
Mj — M, so we can exploit the unitarity relation ^^ |y^^p = 1. Then using 
Oiw = ff|/(47r) and J^i 



\m 



y|2 



((/3)^TOjy, we see that the remaining inte- 
gral is of order Ai^^ where Ai is the larger of the two values M and M\y- 
For the actual numerical evaluation we can also neglect the quark and lepton 
masses. Setting Mw = 80.4 GeV, nitop = 172.9 GeV, ^m„ = 1 eV, and (as an 
example) ((/?) = 400 GeV we get for various values of M the following numbers 



M [GeV] 


/a[10i8 GeV] 


100 


3.1 


150 


2.6 


300 


2.2 


500 


2.1 


700 


2.2 


1000 


2.3 



As we already pointed out in the introduction, these values may be affected 
by higher order QCD corrections, because as is large at small momenta. 

Acknowledgment: We are grateful to Wilfried Buchmiiller for enlightening 
comments on a first version of this paper. H.N. would also like to thank Manfred 
Lindner for discussions. 
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A Weyl-spinor conventions 

Since we heavily use two-spinor notation throughout this paper we here briefly 
summarize our conventions and notations, see [24] for more information. Em- 
ploying the 'mostly minus' metric 77^,^ = diag(-|-l, —1, —1, —1) we define 

a^ = (1 , a^) , a^"" = e'-^e-^a^^. ^ (1 , -a^) 

where e^"^ = £21 = —£12 = — £^^ = 1 and £11 = £22 = 0, hence £Q-y£''''^ — 5^ (the 
definitions for £"'^ are the same). Then we have 













a^'a" + CT'^a^ = 


2 77^" 








as 


well 


as 


the completeness 


relations 












T?!! 








Tra'^CT 


'^ = 277''% 


<. 


^l'-~ 


-'^€5i 






-T U 


irtnei 


'ITU 


Dre 
Tr 


cr^CT^cr^CT-^ = 


2(r,M!^^p^ - 


-77^^ 


^v\ ^ 


77^^77''" - 


i£^^ 


ypX 



where ^'^^^s — i_ Xo relate 2-spinors to 4-spinors we need the Dirac 7-matrices 
and the charge conjugation matrix C, which are given by, respectively, 

A Dirac 4-spinor ^ then decomposes into two Weyl spinors via 



*- :: - .r ^ w^*v-(x",^.) 




where the indices on Weyl spinors are always pulled up and down from the left, 
e.g. (^" = e°'^ipj3 ^ fa = Sa/sf , furthermore (ifia) = 'Pa, etc. The charge 
conjugate spinor is 

so ^ is Majorana if (pa = Xa- Note that in the main text we label the in- 
dependent Weyl spinors by the subscripts L and R, as in (8) and (9), with 
the exception of the neutrino 4-spinor, for which we use different letters i^ and 
N. This is done mainly in order not to encumber the notation with too many 
different letters, although it does not quite conform to standard usage, where 
*L,K are usually defined as the degenerate (projected) 4-component spinors 
i(l =F7^)\['. Finally, we recall that hermitian conjugation inverts the position of 
the (anti-commuting) fermionic operators, i.e. {(px ■ ■ ■ ipY = t/j^ • • • x^V^- 
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